We now use the snp_effects() function... didnt know that existed, but makes sense

I also add grey lines to the amino acid plot of the diplotype clustering functions - it looks a lot neater.

And the haplotype function can now support multiple transcripts.

This is super cool, @sanjaynagi. Thanks very much. I am struggling with bandwidth this week but have asked @leehart and @jonbrenas to take a look so we can merge asap. :)

Awesome, thanks @ahernank

Just adding a corresponding issue for this #819

Many thanks @sanjaynagi - This is great!

Pylint is telling me there's a potential glitch around the following code in the transcript_haplotypes function in the AnophelesHapClustAnalysis class:

        # Get SNP genotype allele counts for the transcript, applying snp_query
        df_eff = (
            self.snp_effects(
                transcript=transcript,
            )
            .query(snp_query)
            .reset_index(drop=True)
        )

Raises problem:

Instance of 'AnophelesHapClustAnalysis' has no 'snp_effects' member

I gather the issue is that AnophelesHapClustAnalysis inherits from AnophelesHapData and AnophelesSnpData, which both inherit from (AnophelesSampleMetadata, AnophelesGenomeFeaturesData, AnophelesGenomeSequenceData) but none of those have snp_effects, which is in AnophelesSnpFrequencyAnalysis instead.

class AnophelesSnpFrequencyAnalysis(AnophelesSnpData, AnophelesFrequencyAnalysis):

Obviously all the tests pass, and the examples in the notebooks appear to work, so I don't quite understand why something like an AttributeError is not already being raised somewhere.

We should probably check to see whether this warning from Pylint is right, anyhow. We might need to modify the inheritance chain, e.g. adding AnophelesSnpFrequencyAnalysis to AnophelesHapClustAnalysis, or move snp_effects into AnophelesSnpData, or maybe there's a different way to call snp_effects, etc.

Also, here is a notebook with a demo of the function https://colab.research.google.com/drive/1HqDJL0cz9gAKaeXdwYjcGXx-X1GYJznL?usp=sharing

Looks like we're still waiting to try using AnophelesSnpFrequencyAnalysis (which has snp_effects) instead of AnophelesFrequencyAnalysis for the AnophelesHapClustAnalysis class, to satisfy Pylint.

I believe AnophelesSnpFrequencyAnalysis will need to be placed before AnophelesSnpData in the list, because AnophelesSnpFrequencyAnalysis inherits from AnophelesSnpData, e.g.

class AnophelesHapClustAnalysis(AnophelesSnpFrequencyAnalysis, AnophelesHapData, AnophelesSnpData):

We might be able to address the potential diamond inheritance problem in a separate issue.

Otherwise, I reckon this looks good to go.